function specres(res)
l = size(res)
k=0;
fr=[res.f0];
jump=[res.jump];
isi = [res.isi];
per = 1000./fr;
l1x=1:1:80;
l1y=l1x;
figure;
plot3(isi, per, jump, '-rs');
grid on;
axis square;
%plot(l1x, l1y);
jlim = find(jump < 30);
jl = [jump(jlim)];
pl = [per(jlim)];
[p, s]=polyfit(pl, jl, 1);
reg = p(1)*l1x+p(2);
%plot(l1x, reg, 'r');

